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Abstract 

We reinvestigated the parameter space allowing a large BR of the 4-body decay of 
the lighter top-squark (fi) accessible at Tevatron Run-II by imposing the lighter 
CP-even Higgs boson mass {ruho) bound from LEP. Important constraints were 
obtained in mSUGRA as well as in the unconstrained supersymmetric models. 
Our results show that the prospect of searching the lighter top-squark via the 
4-body decay mode, in particular the i -\- n — jets + pr signal, is not promising 
in mSUGRA due to the above bound on ruyfi. The existing bounds on 
from Tevatron Run-I and LEP assuming 100% BR of the loop decay of ti are, 
therefore, valid to a good approximation. We also hnd that large BRs of the above 
4-body decay are allowed in the unconstrained model over signihcant regions of 
parameter spaces and the possibility that this decay mode is the main discovery 
channel at Tevatron Run-II is open. We have briefly reviewed the theoretical 
uncertainties in the calculation of ruho and their consequences for the constraints 
obtained by us. We have commented upon, with illustrative examples, how the 
above parameter space is affected if future experiments push the Higgs boson 
mass bound upward. 
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1 Introduction 


The mass of the lighter top-squark (ti), superpartner of the top-quark in the standard model 
(SM) may turn out to be below the top-quark mass. This happens in a wide region of the Su¬ 
persymmetry (SUSY) parameter space both in the minimal Supergravity (mSUGRA) model 
and in the Minimal Supersymmetric Standard Model (MSSM). It could even be the next- 
to-lightest supersymmetric particle (NLSP), the lightest neutralino (y)*) being the lightest 
supersymmetric particle (LSP) by assumption. Being the lightest among the strongly in¬ 
teracting sparticles with a relatively large cross-section compared to other sparticles, the 
lighter top-squark may show up in direct searches at Tevatron Run-II. In our convention, 
ti = R cos 9j + tR sin 9j and t 2 =-tL sin 9j+t^ cos 9j, where 9j is the mixing angle and G is 
the heavier top-squark. The top-quark mass {rrit) appears in the off-diagonal elements of 
the top-squark mass matrix. So, the physical top-squark masses and the mixing 

angle depend on the input value of the top-quark mass. 

In the lighter top-squark-NLSP (ti-NLSP) scenario {i.e. , m~o < < the mass of any 

other sparticle), U has few allowed decay modes in the R-parity conserving (RPC) SUSY 
model. They are the Flavor Changing Neutral Current (FCNC) induced loop decay, U —> cXi 
[1], the 4-body decay, ti bxiff', where //' = u,d,c, SjijUe [2], the tree-level two body 
decay, ti —> txi and the three body decay, ti bWxv If ^ nit + production 

cross-section of ti is rather low for the Tevatron experiments. The ti bWxi mode may be 
suppressed if the coupling WxiXi is small or < Mw + rn-o +mb. This may happen if Xi 
is bino like as in the mSUCRA model. However, in this report we have considered only those 
parameter spaces where the last two decay modes are either kinematically or dynamically 
suppressed. 

So far most of the ti-NLSP searches at Tevatron Run-I and at the Large Electron Positron 
(LEP) collider have been performed assuming 100% BR of the loop decay^, leading to the 
signal: jets -|- missing energy {^). The most stringent constraint comes from the Tevatron 
Run-I experiments [5]. However, the above assumption becomes unrealistic, if ti has other 
decay modes with comparable widths. This possibility exists in a large region of the SUSY 
parameter spaces in both the models mentioned above [2, 6]. The 4-body decay and the 
loop decay in RPC models may have decay widths of the same order of magnitude in a wide 
region of the parameter space [2]. 

An interesting feature of all supersymmetric models is the prediction of the existence of 
at least one light Higgs boson [7, 8]. So far the LEP and Tevatron Run-I experiments found 
no signihcant evidence in favor of a light Higgs boson. The direct searches from the four 
LEP experiments concluded that the lighter Higgs boson must be heavier than 114.4 CeV at 
95% C.L. [9]. Generally this bound depends on the value of sin^(/5 — a), where f3 and a are 
the ratio of the vacuum expectation values (VEVs) of the two Higgs doublets and the mixing 
angle in the CP-even Higgs boson mass matrix respectively. However, the Higgs boson mass 

^Recently the D0 collaboration [3] and the ALEPH collaboration [4] looked for the ti b^ff' channel 
with special assumptions about its BR. 
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bound is 114.4 GeV when sin^(/3 — a) ~ 1. We have found that this factor is indeed ~ 1 for 
the parameter space interesting for us. This Higgs boson mass bound is directly applicable 
in SUSY when the CP-odd Higgs boson {Ma) mass is very large compared to Z-boson mass 
{Mz), i.e. , Ma » Mz, commonly known as the decoupling limit. 

In SUSY the mass of the light Higgs boson at tree-level is given by, 


mlo 


Ml + Ml- ^j{Ml + Miy - 4M|Ml cos^ 2(3 


( 1 ) 


Eqn.l shows that rriho ~ Mz at the tree-level. But radiative corrections involving top 
quark and squarks (t) in the loop indicate that m^o grows like CpNcCrnl where Gp \s the 
Fermi coupling constant, is the color factor and C is a model dependent loop factor. 
It is clear that the relatively large experimental error in top mass may lead to a sizable 
uncertainty in the Higgs boson mass prediction. The loop corrections from bottom quark 
and squarks (6) are also sizable for large values of tan/? and /i [10], where /i is the Higgs 
mass parameter in the superpotential. We shall discuss the Higgs boson mass as a function 
of the SUSY model parameters in section 2. 

Theoretically the mass of the lighter CP-even Higgs boson is bounded from above: rriho ~ 
135 GeV. This bound has been obtained by including radiative corrections up to two-loop 
level [8, 11, 12, 13, 14, 15]. So, rriho is expected to he in the range 114.4 GeV ~ ~ 135 GeV. 

Due to the yet unknown higher-order corrections the theoretical error on the lighter CP-even 
Higgs boson mass consists of several pieces [16, 17, 18, 19, 20]. They are the momentum- 
independent two-loop corrections, the momentum-dependent two-loop corrections, the higher 
loop corrections from t/t sector etc. Taking into account of all shorts of uncertainties the 
intrinsic error has been estimated to be 6mho ~ 3 GeV However the magnitude of this 
uncertainty crucially depends on the SUSY parameter space. In view of the above theoretical 
errors the viability of a parameter space will be judged in this work by using the following 
two constraints on the calculated Higgs mass: (i) niho ~ 111.4 GeV (the weaker bound) 
{ii) ruho ~ 114.4 GeV (the experimental bound). The first bound obviously leads to more 
conservative results as will be shown in section2. 

We now make a few comments on the ti decay widths considered in this paper. In the 
approximation of neglecting the charm qnark mass, cr does not mix with top-squarks . The 
mixing of cp with the top-squarks mass eigenstates result in a redehned lighter top-sqnark 
given approximately by ii=ii + ecp, where 

ALCOs6<y-Ai?sin6'y 
e 2 2 ■ 

mi — mi^ 

t\ 

The parameters A^^/? are given by 


Ar, = 


a 
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^Inclusion of the full two-loop corrections as well as the leading higher loop corrections is expected to 
yield a reduced uncertainty Srriho ~ 0.5 GeV [16, 18] 
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( 4 ) 


where mg^, rnin^ ^b are the left-handed scharm, right-handed sbottom, down- 

type Higgs scalar and trilinear sbottom mass parameters respectively. Mqut is the Grand 
Unification theory (GUT) scale in SUSY model. 

Assuming proper electro-weak symmetry breaking (EWSB), the Higgs scalar mass 
squared parameter (in Eqn.3) can be written as 



( 5 ) 

( 6 ) 


r at , AT \ 2 2 ^ -swNn + cwNu] ^ ^ 

/li — V 2 -{cwNii -h swNu) + o ~ Q -Sty 1 - (7) 

.o \/ O / Cw-Sw 

Here Yu and Y 12 denote the elements of the neutralino mixing matrix projecting the LSP on 
to the photino {j = 1) and zino {j = 2) states. After renormalization one obtains the residual 
logarithmic terms in Eqns.3 and 4 which lead to a large contribution in the amplitude. 
Adding various diagrams contributing to the loop, a divergent term is obtained in general. 
In principle this divergence must be cancelled by adding appropriate counterterms. The 
dominant effects of renormalization are expected to come from the logarithmic divergences 
caused by soft breaking terms at Mqut in Supergravity (SUGRA) models. Thus in practice 
a large logarithmic factor ln(MQUT/''^w) ~ 65 is included in the decay amplitude. 

The 4-body decay mode U — bxiff', proceeds through various diagrams [2]. An approx¬ 
imate analytical expression for the 4-body decay amplitude can be found in [2]. However, 
the decay rate has been been calculated in [2] without any approximation by taking into 
account all diagrams and the interferences among them. We thank the authors of [2] for 
providing their FORTRAN code. The packages SDECAY [21] and CalcHEP v2.1 [22] are now 
available for calculating the decay width in the context of MSSM. We have cross checked the 
FORTRAN code with CalcHEP v2.1. 

In general, in any diagram the virtuality of the exchanged particle in the propagator is 
important for getting an appreciable contribution. This gives an idea of the parameter space 
where the 4-body width is signihcant. It has been checked that for the ti mass range and the 

■^For SUSY breaking mechanisms other than gravity mediation, such as gauge mediation, the SUSY 
breaking scale can be much lower than the GUT scale, which lead to a somewhat smaller value of F in 
Eqn.6. 
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model parameters used in this analysis the chargino and/or the slepton mediated diagrams 
are the dominant ones and the 4-body width is indeed comparable to the loop decay width. 

The 4-body and the loop decay of the ti-NLSP are of the same order in perturbation 
theory, i.e. O(a^). So, they may compete with each other. Our prime objective is to hgure 
out the SUSY parameter space where the 4-body BR may compete with or overwhelm the 
loop decay BR. In this perspective a few more points are noteworthy. 

• As discussed above the large logarithmic factor in the width (see, Eqns.3 and 4) pre¬ 
supposes the existence of a grand unihcation scale. Moreover, it is also assumed that 
MSSM is the correct theory all the way up to Mgut- If this is not the case this factor 
can be rather small. This would enhance the 4-body BR compared to the typical values 
presented in this paper. 

• If the ti is purely right-handed, the amplitude involves only Ar in Eqn.2. This will 
suppress the amplitude for small values of the trilinear coupling Ai, (in Eqn.4). Addi¬ 
tional suppression may come if the (common) SUSY-breaking scalar mass of the hrst 
two generations is large (Eqn.2). 

• Any scenario with tan^j ~ Al/Aji immediately leads to a negligible e (Eqn.2) . 

• From Eqns.3 and 4 we note that the mixing between ti and cl is proportional to 
tan^/5. It is, therefore, relatively suppressed at small tan/? which is favourable for a 
large 4-body BR [2, 6]. 

It has already been pointed out that the ti searches in the 4-body leptonic decay channel 
is especially promising for 3 ~ tan/? ~ 5 at the upgraded Tevatron [6]. However, it is also 
well known that the radiatively corrected rriho increases with tan j3. So, what happens to the 
parameter spaces favourable for ti searches in the 4-body decay channel if the Higgs boson 
mass bound is invoked? We are addressing this issue in this paper. 

The input value of top-quark mass in this paper is important for two reasons. Firstly, 
it changes the parameter space where the ti-NLSP condition is satished. Consequently the 
4-body decay BR will be affected as well for a particular choice of SUSY model parameters. 
Secondly, as discussed above the, radiative corrections to the Higgs boson mass is very 
sensitive to m*. 

Before delving into the numerical discussions let us briefly review the current status of 
the top mass measurement. From the published Tevatron Run-I and preliminary Tevatron 
Run-H data, the present world average for the top-quark mass is rrit = 172.7 ± 2.9 GeV [23]. 
Since there is considerable uncertainty in mt we shall consider in this paper two representative 
values: (i) 771^=172.7 GeV which is the central value of the Run-I and Run-H average and {ii) 
mt=178.0 GeV which is on the higher side of the currently allowed range ( approximately 
2a away from the present central value ). The second choice predicts relatively large rrtho 
for hxed values of other SUSY parameters and hence, leads to weaker constraints on the 
parameter space allowed by the bound on niho. 
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The plan of the paper is as follows. In section 2, we begin with a discussion of the nu¬ 
merical tools used for calculating the sparticle masses and the lighter CP-even Higgs boson 
mass. The sub-sections 2.1 and 2.2 are devoted to the phenomenological discussions of the 
lighter top-squark 4-body decay modes respectively in mSUGRA and MSSM models consis¬ 
tent with the lighter CP-even Higgs boson mass limit from LEP experiments. Consequences 
of the theoretical uncertainties in the rriho prediction are also analyzed. We shall conclude 
in section 3. 


2 Numerical Analysis : 

We begin with representative parameter spaces in which the fi-NLSP condition is realized. 
We then draw the contours of its total 4-body decay BR in this parameter space. Finally 
we impose the lighter CP-even Higgs boson mass constraint on this parameter space. The 
sparticle spectrum in mSUGRA has been calculated using ISAJET v7.48 [24]. The sparticle 
spectrum in MSSM has been calculated using our own codes. The Higgs boson mass has 
been calculated using SUSPECT v2.33 [25] in both the models. This code uses the full one- 
loop and dominant two-loop contributions to the self energy corrections in the SUSY Higgs 
boson mass matrix in the Dimensional Reduction {DR) scheme 

The supersymmetric radiative corrections to the particle (top, bottom, tan) and spar- 
ticles masses (squarks and gauginos) have also been incorporated. We have set the follow¬ 
ing input values in SUSPECT v2.33 at the weak scale {Mweak)- ^ 27 ^ 934 ’ «s=0.1172, 

mt(pole)=172.7 or 178.0 GeV, mT-(pole)=1.777 GeV, mfe(mfe)=4.25 GeV, sing^=0.2221. 

The Ma and /i are treated as inputs in MSSM in this paper. All masses and mass 
parameters in this paper are in GeV. 

The magnitude of /i can be constrained from the radiative electro-weak symmetry break¬ 
ing (REWSB) conditions which produce the observed Z-boson mass. Minimizing the Higgs 
potential with radiative corrections one obtains 


and 


-Mi = 




tan^ /5 


tan^ f3 


/i' + A 


R 


fiB = 


m 


Hi 


+ 


-I- 2/i^) sin 2/5 . 


( 8 ) 

(9) 


, where are the soft Higgs mass parameters at M^^ak and A/j is the radiative 

correction which is given in [24] and [25] . This REWSB criterion restrict the \^B\ and 
\fj,\ at Mweak and leaves the sign{fi) as a free parameter in mSUGRA. One can trade the 

^Other public codes, for example, FeynHiggs [26] and CPsuperH [27] are also available for calculating the 
properties of the Higgs sector in various SUSY models. 
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parameter B for tan/? (ratio of Higgs VEVs). Eqns.8 and 9 are also valid in MSSM bnt for 
arbitrary soft masses. 

The soft SUSY breaking terms cannot be arbitrary at M^^ak in mSUGRA. The soft terms 
are renormalization group equation (RGE) driven from high scale inputs like the common 
scalar mass, see sub-section 2.1. The trilinear couplings at My^^ak can also be computed from 
a common input value at Mquj' using the RGE. The soft terms are treated as independent 
inputs in MSSM which we shall discuss in the sub-section 2.2 . 

However, in both the models one can show that a deep charge-color-breaking (CGB) 
minima [28, 29, 30] appears in the scalar potential at My,eak unless 


I Ap ^ 3(m2- + m}^ + m]j^ + /i^) 

^ 3(m^-^ + + m]j^ + (10) 

~ 3(m^-^ + ml-^ + + f/) 

where At, and A^ are the trilinear soft parameters of the top-squark, bottom-squark and 
tau-slepton sector respectively. 

In our analysis we have considered the ranges of model parameters such that the CGB 
minima [28, 29, 30] of the potential do not arise. Moreover all sparticle masses are required 
to satisfy the recent experimental bounds. We now make a few comments on the obtained 
ti mass limits. 

The published mass limits on ti not only depend on the production cross-sections, they 
also depend on AM = — misp, where LSP can be either Xi oi' In addition these 

limits depend on the BRs of various decay channels kinematically allowed. Moreover, the 
limits from LEP also depend on Oj. 

Most of the mass limits from Tevatron [31, 32, 33, 34] and LEP [35, 36] have been 
obtained by assuming 100% BR of the ti cXi or the ti bxt decay. However, the 
presence of other competent decay mode, for example the 4-body decay mode emphasized in 
this paper, substantially changes these limits. The CDF limit [32] with the above assumption 
has been included in the analysis of section 2.1. 

Recently the D0 collaboration [3] and the ALEPH collaboration [4] put limits on 
assuming complete dominance of the channel ti bx^lvi and coexistence of the loop decay 
and the 4-body decay respectively. From Fig.4b of [4] the absolute lower limits of is 
found to be 63 at 95% C.L. . This limit is obtained for AM = 5, BR{ti —> cx5’)=22 % and 
BR(ti —> b^lvi) =55 % and 6*y= 56°. The my limit becomes stronger when BR{ti —> cXi) 
is increased from 22%. It may be noted that the assumed BRs in the above works are not 
realistic. For example it was shown in [6] that the BR of the leptonic {i = e and fi) 4-body 
decay mode is unlikely to exceed 20% since competing hadronic 4-body decay modes are 
inevitably present. In a different approach [37] model independent limits on the product of 
branching ratios BR{ti be'^Ve^) x BR{t\ bqq'Xi), where q and q' correspond quarks of 
all flavours kinematically allowed, were obtained as a function of my using Tevatron Run-I 
data. 
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The LEP-II collaboration puts the limit ~ 96 [38] assuming complete dominance of 

hi bxi decay. The following constraints on sparticle masses have also been included in the 
analysis of this paper: m~ ~ 86, m~ ~ 43.7, m~^ ~ 94.9, m~^ 7 99.4 [39] and m~± ~ 103 [40]. 

2.1 The mSUGRA model: 

The input parameters of mSUGRA at Mgut are the common scalar mass (mo), the common 
gaugino mass (mi/ 2 ), the trilinear scalar coupling term (Aq); tan and sign of /i (the higgsino 
mass parameter) [41]. The magnitude of /i is hxed by the REWSB condition [28, 42] and 
depends on the input top quark mass®. 

Since the NLSP nature of the ti is very crucial for our work, we have demarcated different 
regions in the Rq ~tan (3 parameter space where ti becomes the NLSP. A detailed discussion 
is given in sec.II of [6] (see Fig.l and 2 ^ of [6]), where it was assumed that mt=175.0 . Similar 
information is given in Fig.l and Fig.2 in this paper for mt=172.7 and mt=178.0 respectively. 
Several contours of the combined BR of all 4-body decay modes are also shown. In order 
to identify the parameter space allowed by the m^o constraints, we present the contours 
of rriho ~111.4 (the weaker bound) in both the hgures and m/jo ~114.4 (the experimental 
bound) in the right panel of Fig.2 only. 

The choice of other mSUGRA input parameters used are given in the hgure captions. 
These choices are guided by the following facts. If uiq » mi/ 2 , the gauginos are naturally 
light. This disfavours the ti-NLSP scenario. For mi /2 >> mo, on the other hand, one 
may require large |Ao| for fhis scenario. However, such high values of |Ao| may violate of 
the GGB conditions, see Eqn.lO. So, in our analysis we focused on uiq ~ mi /2 regions for 
mSUGRA. However, the ti-NLSP is quite common in larger regions in the parameter space 
of the unconstrained MSSM since the soft breaking parameters are arbitrary. This model 
will be analyzed in the next section. 

For each tan/3 there is a range of |Ao| ( l^^oLm < l^^ol < l^olmax ) ^^gut which 
corresponds to the fi-NLSP (see Fig.l and 2). Otherwise fi or y® m xf happens to be the 
NLSP and U decays into other 2-body or 3-body decay channels. 

For |Ao| > \7^o\max^ becomes the LSP which is forbidden. On the other hand if 
|Ao| < 1^01 min fo^ given set of mSUGRA parameters, ti becomes heavier than the lighter 
chargino and the ti-NLSP condition is violated. We, therefore, concentrate on |Ao| ~ I ^ 0 1 min 
so that the chargino virtuality is small [2, 6] and consequently the 4-body BR is appreciable. 
The off-diagonal terms in the top-squark mass matrix have the form [At -:j;^^). So as long 
as At at Mweak and fi have opposite signs the NLSP scenario is favored. 

For mi=172.7, the ti-NLSP is realized in a large region of the parameter space (see 
Fig.l). As expected the total BR of the 4-body decays is appreciable for relatively low value 
of tan/3 [2, 6]. The region above the rriho ~ 111.4 contour cannot be excluded due to the 

®The top-quark mass is very crucial for calculating ttiHu and Aji in Eqns.8 and 9. 

^One can easily find m~o and m~± from the parameter set used in the figures. 

Xi Xi 
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theoretical uncertainties mentioned in the introduction. It should be noted that even with 
the weaker bound on rriho the total 4-body decay BR is ~ 5% for mo = 200, mi /2 = 145 (see 
the left panel of Fig.l). For mo = 140, mi /2 = 180, on the other hand, the total 4-body decay 
BR can be as large as 40% (see the right panel of Fig.l). However, the whole parameter 
space in Fig.l is disfavoured if the experimental bound is imposed. 

A similar analysis for mt=178.0 is shown in Fig.2. In the left panel of Fig.2 with mo = 
200, mi /2 = 145 we have some allowed region if the weaker bound on mhO is imposed. Still 
the total 4-body decay BR can be at most 20%. In the right panel of Fig.2 with mo = 140, 
mi /2 = 180 we have some allowed region even if the experimental bound is imposed. However, 
the total 4-body decay BR can be at most 20%. Of course with the weaker bound larger 
4-body decay BRs ( ~ 60%) are allowed. It should also be noted that for a hxed mo, mi /2 
cannot be increased arbitrarily because that would lead to other instabilities of the scalar 
potential [28, 29, 30, 43]. 

We have shown the variation of mj^{mho) as a function of Ao in the left (right) panel of 
Fig.3 for mt=178.0 for three values of tan 13. Other model parameters are given in the hgure 
caption. The ti-NLSP criterion is not imposed in this particular hgure. However, the range 
of mj^ satisfying the ti-NLSP condition can be easily found out by comparing with the right 
panel of Fig.2. Then one can read off the lighter chargino and neutralino masses from the 
left panel of Fig.3. A ti having mass in the range shown in Fig.3 may be copiously produced 
at Tevatron Run-H. From the right panel it is clear that tan /5=5 is marginally allowed only 
if the weaker bound on m^o is required. On the other hand if tan (3^9 then m^o is consistent 
even with the experimental bound, but the total 4-body decay BR is negligible as can be 
seen from the right panel of Fig.2. 

The scatter plots in Fig.4 demarcate the ti-NLSP region consistent with the bound on 
mho in the mo — mi /2 parameter space for positive p, mt=178.0, tan/5=5 (left panel) and 
tan/5=9 (right panel). In the left panel of Fig.4 with Ao="630 we have some region allowed 
by the weaker bound on m^o. The total 4-body decay BR can be as large as 50%. In the right 
panel of Fig.4 for tan (3=9 we have some allowed region even after imposing the experimental 
bound on m/jo but the total 4-body decay BR can be at most 10%. With the weaker bound 
on m/jo the whole ti-NLSP region is allowed. In the regions with m^o — 111.4 the total 
4-body BR is <5%. 

We have repeated the analysis for mt=172.7 using the same parameters as in Fig.4. In 
this case the entire ti-NLSP region in the mo — mi /2 parameter space is disfavoured for 
tan/d=5 even for the weaker bound on mho ~ 111.4. For tan/5=9 on the other hand almost 
the entire parameter space is allowed by the weaker bound on mho. However, the total 
4-body BRs can be at most 10% for tan (3=9. 

In [6], we have studied the H search prospects at Tevatron Run-H using the one lepton 
with 2 or more jets accompanied by a large amount of missing energy signal assuming 
mt=l75. This signal topology arises from the pair production of the lighter top squark 
followed by the leptonic 4-body decay of one lighter top-squark (ti bxiiue) while the other 
decays hadronically (ti bx^qq'). We have shown that, BR(ti bx^iui) ~ 20%, where 
i = e and p (see Figs.3 and 4 of [6] ), is adequate for detection at Tevatron Run-H. From 



the present analysis we find that such large leptonic BRs cannot be realized at least in the 
mSUGRA model for mt=172.7 and not even with the rather conservative choice mt=178.0 
irrespective of the rriho bound invoked. Our analysis also show that the ti mass limits from 
LEP and Tevatron Run-I assuming 100% BR of the loop decay of ti are approximately valid 
in the context of mSUGRA if rrit is close to its current central value (172.7) irrespective of 
the uncertainties in the predicted m/jO. However, in view of the uncertainties in rrit and in 
the calculated m^o one cannot conclusively exclude the possibility of a relatively large total 
BR of all 4-body decay modes (see, e.g. , the right panel of Fig.2). 

We have shown some points in the parameter space excluded by the GDF limits on 
in Fig.l (left panel). Fig.2 (left panel) and Fig.4 using the heavy dotted squares. For each 
of these points m~o — 51, ~ 102 and 4-body BR 5 10%. It is heartening to note that 

even the modest statistics from Tevatron Run-I can already put constraints on Aq which is 
rather difficult to constrain from other measurements. However, the GDF excluded points 
are already ruled even by the weaker bound rriho ~ 111.4 in the context of mSUGRA. One 
notable exception is found in the right panel of Fig.4. Here we find that a small region of 
parameter space allowed by the weaker bound on m/jO is ruled out by the GDF constraint. 


2.2 The MSSM model: 


In the MSSM with a universal gaugino mass, the SU(2) gaugino mass parameter (M 2 ), /i and 
tan f3 completely describe the neutralino and chargino masses and mixings. Depending on the 
relative magnitudes of M 2 and /i, the model is categorized in three different scenarios. They 
are the higgsino (M 2 >> /i), the gaugino (M 2 « jj) and the mixed (M 2 ~ /i)scenarios. The 
lighter chargino (xi^) and the two lighter neutralinos (x5 and y®) have approximately the 
same mass ( ~ /i) in the higgsino scenario. Thus it is difficult to accommodate the ti-NLSP 
without fine adjustments of the parameters. The ti happens to be the NLSP in large regions 
of the parameter space in the gaugino and the mixed scenarios. In this paper we studied 
the mixed scenario. However, the main results and conclusions are qualitatively valid in the 
gaugino scenario. 

To calculate and rrihO in MSSM, we use the parameter (instead of Aq) as input, 
where 


^ _ {At - fj, cot P) 


( 11 ) 


rriuL and are the soft SUSY breaking mass parameters for the first generation left and 
right-squark respectively at the weak scale. 

In our analysis the first generation left and right-squark soft masses are assumed to 
be the same, i.e. , m~^ = m~^ = mj = m~. The charged slepton sector is also treated 
similarly, i.e. , m~^ = m~^ = mj, but the universality of squark and slepton masses is not 
assumed. These mass patterns can be realized in some variations of the mSUGRA model. 

For the third generation squarks the above equality may not hold even in mSUGRA due 
to, for example, RGE driven by the Yukawa terms. However, the SU (2)x, invariance require 
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and neglecting the D-term contribntions. In this paper the right and 

left-sqnarks soft masses of the third generation are assnmed to be eqnal to m~. The physical 
masses are calcnlated after incorporating the corresponding fermion masses, the D-terms etc 
in the mass matrix. 

For a given fi, tan (3, m~ and Xr, the trilinear soft breaking mass parameters {At) of the 
top-sqnark sector can be calcnlated from Eqn.ll. The same parameter for bottom-sqnark 
{Ab) and tan-slepton {Ar) sectors are treated as free inpnt parameters. Using these inpnts 
the third generation sqnark and slepton masses and mixing angles can be calcnlated. 

We have nsed the CP-odd Higgs boson mass (M^) and /i as the free parameters in 
SUSPECT v2.33 to calcnlate the lighter CP-even Higgs boson mass. 

The off-diagonal terms in the top-sqnark mass matrix have the form mtXrfn-. So, the 
masses and the mixing angle of the top-sqnark sector are controlled by these parameters. For 
negative Xr, the and are exactly the same as that for positive Xr. The mixing angle 
depend on the relative magnitnde of the diagonal and off-diagonal terms in the top-sqnark 
mass matrix. 

The my (m/jo) as a fnnction of Xr has been shown for mt=172.7 in the left (right) panel 
of Fig.5 for fonr valnes of tan/? (see the hgnre caption for other inpnt MSSM parameters), 
my is insensitive to m* and tan (3. 

The lighter CP-even Higgs boson mass has a periodic dependence on Xr. For a given 
Xr, if tan/? increases the Higgs boson mass also increases. If tan/5=5 (15) and mi=172.7, 
rriho ~ 114.4 for Xr ~2.0 (1.5) and m/jO ~ 111.4 for Xr — 1.75 (1.30). It is clear from Fig.5 
that both the ti-NLSP condition and the experimental bonnd on rriho may be satisfied for 
large positive valnes of Xr. In onr snbseqnent analyses we shall only consider snch valnes of 

Dne to the large nnmber of free parameters in MSSM we varied two important parameters 
at a time keeping the rest hxed. 

• Xr — tan (3 

We have demarcated the ti-NLSP regions in Xr — tan [3 plane in the left (right) panel of 
Fig.6 for mt=172.7 (178.0). Several contonrs for hxed valnes of the total HR of 4-body decays 
are also shown. Comparing with Fig.5 we hnd that the mass of the fi-NLSP is kinematically 
accessible at Tevatron Rnn-H. For each tan f3 there is a range of Xr ( < \Xr\ < 

) satisfying the fi-NLSP condition ( recall the Aq ranges in Figs 1 and 2 ). 

From each HR contonr it follows that when Xr is increasing smaller valnes of tan (3 are 
reqnired. For large Xr, my decreases and the virtnality of the exchanged chargino and/or 
slepton becomes large and tends to snppress the 4-body decay width. In order to keep the 
HR hxed tan (3 mnst be correspondingly rednced. 

We have also plotted several contours for m/jo ~ 114.4 and a few other Higgs boson 
masses. Comparing the left panel of Fig.6 with the right panel of Fig.5 we hnd that for 
m/jO ~ 111.4 the ti-NLSP condition is violated. 

As has already been noted, for mt=178.0 the constraints imposed on the parameter space 
by the bound on m/jO are less severe (see the right panel of Fig.6). The total 4-body HR can 
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be as large as 90% for tan /9 ~ 9 for both choices of rrit in Fig.6. Even a t tan ~ 24 the total 
4-body BR is small but not negligible (~ 10%). In this case it will be interesting to search for 
hadronic 4-body decays of ti at the Large Hadron Collider (LHC) or the International Linear 
Collider (ILC). Even if tan/? ~3, rriho ~114.4 is allowed. Consequently almost the entire 
parameter space in the right panel of Fig.6 is allowed. Such low values of tan/? naturally 
guarantees large total BR of the 4-body decay modes. Moreover we hnd that the BR of 
leptonic 4-body decays can be as large as 20%, required for an observable signal at Tevatron 
Run-H [6]. 

It is expected that the upcoming collider searches will either discover the Higgs boson 
or push its mass bound upward. Moreover, the theoretical uncertainties on the Higgs boson 
mass is likely to be signihcantly reduced {Sniho ~ 0.5) [16, 18]. We would now like to comment 
on the impact of such possible improvements in our understanding of m^o on ti searches. 
For illustration we have plotted contours of niho ~ 116.0, 118.0 and 119.0 (118.0, 120.0 and 
122.0) for mt=l72.7 (178.0) in the left (right) panel of Fig.6. The total 4-body BR ~ 30% 
can be completely (nearly) ruled out by an improved Higgs boson mass bound, rriho ~ 119.0 
(122.0) if mt=l72.7 (178.0) for our choice of other SUSY parameters. 

• Xr — Ma 

We have shown the ti-NLSP regions, a few 4-body BR and rriho contours in the Xr — Ma 
plane for mt=172.7 and tan/5=5 (9) in the left (right) panel of Fig.7. The other MSSM 
model parameters used are mentioned in the hgure caption. 

The lighter top-squark mass is almost independent of Ma if the other MSSM parameters 
are hxed. So, the H-NLSP region remains the same when Ma is varied for hxed Xr. The 
input value of Ma can affect the CCB condition (see Eqns.5, 8 and 10) through the soft 
breaking Higgs boson masses. This disallows some regions of the parameter space for large 
Xr- However, the CP-even Higgs boson mass depends on Ma- On the other hand Ma has 
a nontrivial impact on the loop decay width, see, Eqns.3 and 5 or Ref.[l]. For a hxed Xr 
if Ma increases the loop decay width also increases. However, the 4-body decay width is 
insensitive to Ma- So, for large Ma the total 4-body BR decreases for a given Xr. This is 
very transparent in Fig. 7 for relatively large values of Xr corresponding to large virtuality 
of the exchanged particles in the 4-body decay amplitude. 

We hnd that the total 4-body BR (leptonic BR) can be as large as 90% (20%) in signif¬ 
icant regions of the parameter space for tan/5=5. If the experimental bound on the Higgs 
boson mass is invoked then the total 4-body BR (leptonic BR) ~ 90% (20%) is allowed if 
Ma ~ 300 (400) and Xr 5 2.17 (2.08) for tan (3=h (see the left panel of Fig.7). We have shown 
TTihO ~ 111.4 contour in the left panel of Fig.7. It is clear that almost entire parameter space 
is allowed by the weaker bound on the Higgs boson mass. It has also been checked that the 
whole ti-NLSP region is disfavoured for m^o -117.4. 

For tan/3=9, in the entire fi-NLSP region in Fig.7 corresponds to rrihO ~ 114.4. However, 
the total 4-body BR (leptonic BR) ~ 90% (20%) occurs only in restricted regions, i.e. , 
Ma ~ 500 and Xr 5 2.025 as can be seen from the right panel of Fig.7. 

Since the lower bound on Higgs boson mass may improve in future, we have shown 
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several rriho contours, e.g. , 115.0 and 115.7 (116.0, 117.0, 117.4, 117.8, 118.4 and 119.0 ) for 
tan/5=5 (9) in the left (right) panel of Fig.7 to illustrate the consequences for the BR under 
study. 


3 Conclusions 


The lighter top-squark turns out to be the NLSP in a signihcant region of the supersymmetric 
parameter space both in constrained and unconstrained models like mSUGRA and MSSM 
respectively. The fi-NLSP has two competing decay modes for the top-squark masses within 
the reach of the Tevatron Run-II. They are the loop induced ti cXi and the 4-body decay, 
ti bxiff- These two decay modes may have comparable decay widths over large regions 
of parameter spaces in both the models [2, 6]. However, for low tan/? the loop decay width 
is suppressed which in turn enhances the 4-body decay BRs. On the other hand the current 
bound rriho ~ 114.4 from LEP disfavours low values of tan/?. It is therefore important to 
identify the parameter spaces allowing large total BRs of the 4-body decays after imposing 
the rriho constraint. 

The theoretical uncertainties in the calculated ruho using any currently available code is 
~ 3 due to the yet unknown higher-order corrections. In view of this we have assumed that 
realistically the bound on the calculated ruho can be relaxed to rriho 7 111.4. 

The masses of ti and rriho crucially depend on the input value of m*. We have considered 
the current world average central value of rrit which is mt=172.7 as well as mt=178.0, which 
is approximately 2a away from the central value and yields somewhat weaker constraints. 

We have demarcated the ti-NLSP regions in the mSUGRA model for mt=172.7 and 178.0 
in the Aq — tan/? ( Figs.l and 2) and mo — mi /2 planes (Fig.4). We found that for a given 
tan/5, a limited range of Aq with Aq ~ |Ao|mm allows large 4-body BRs. For demonstration 
we have presented in the Aq — tan (3 plane a few 4-body BR contours for two representative 
choices of uiq and mi /2 (see sub-section 2.1 and the hgure captions). For mt= 172.7 all 
parameter spaces with the total BR of the 4-body decays > 10 (40)% are disfavoured even 
if the weaker Higgs boson mass bound is invoked (see the left (right) panel of Fig.l). The 
experimental bound on m/jo rules out any parameter space with numerically significant BR 
of the above decays. For mt= 178.0 the experimental bound on rriho leaves no room for total 
BR> 10 (20)% (see the left (right) panel of Fig.2). The weaker bound, however, allows BRs 
as large as 30 - 50 % in a small region of the parameter space (the right panel of the above 
hgure). 

Gomplementary information in the mo — mi /2 plane, for Ao=-630, sign(/i) > 0 and 
mt=178.0 are presented in Fig.4. For tan/5=5 the weaker constraint on m/jo allows BRs as 
large as 60% (see the left panel) in a small parameter space, while the experimental bound 
on m/jO practically rules out the possibility of a sizable BR of the decay channels of interest. 
For tan j3 = 9, the total 4-body BR can be at most 10% due to the experimental bound on 
Higgs boson mass. 
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We have also excluded small regions in Aq — tan/? and mo — mi /2 plane from the CDF 
limits on the mass of ti in the jets+^ channel. In most cases, however, these constraints 
are superseeded even by the weaker bound on rriho. 

In conclusion we note that in mSUGRA the uncertainties in m/jo and mt leave open the 
possibility that the total BR of the 4-body decay of the ti-NLSP is sizable although the 
parameter space corresponding to this scenario is severely squeezed by the bound on m/jO. 
Perhaps the main future interest in this decay channel in the context of mSUGRA would be 
to observe it as a rare decay mode at the LHC or the ILC. Moreover, the existing bounds 
on obtained by assuming 100% BR of the loop decay are approximately valid in most of 
the parameter space. 

In MSSM we have used the parameter X,.(see Eqn.ll) as an input. The interplay of the 
ti-NLSP condition and the Higgs boson mass constraints can be easily followed by varying Xr 
(see Fig.5). One of the important conclusions is that sizable total BR of the 4-body decays 
of ti can occur for the stronger constraints rrit = 172.7, m/jo ~ 114.4 and intermediate values 
of tan (3 (10 - 15) (see Fig.6). For example, this BR can be as large as 90% in the X^. — tan (3 
plane even for tan/? ~ 9, rrihO ~ 114.4 and m* = 172.7 and 178.0. For lower values of tan/? 
the leptonic 4-body BR ~ 20%, which can lead to an observable signal at Tevatron Run-H 
[6], is also consistent with the above constraints. The total BR can be restricted to the level 
of 10 - 20 % only if the bound on m/jo is signihcantly improved and rrit is determined with 
better precision. Complementary information leading to similar conclusions are presented 
also in the Xy. — Ma plane (see Fig.7). 

In conclusion we emphasize that the bounds on from Tevatron obtained by assuming 
the complete dominance of the loop decay of ti requires revision if we want these bounds 
to be model independent. This is true even though the strong constraint on ruho from LEP 
squeezes the parameter space which favours the competing 4-body decays of ti. Based on 
our earlier analyses we conclude that the 4-body mode could still be the main discovery 
channel for ti at both Tevatron Run-H and LHC. As future strategies for the discovering ti 
are planned, attention should be paid to the lighter CP-even Higgs boson search. If its mass 
bound is signihcantly improved, the interest in the 4-body decay would be to look for it as 
a rare decay. It is expected that the Higgs boson may show up at an early stage of LHC the 
experiments and, at least, some preliminary information on its mass will be available. This 
information will hnally decide the viability of ti search in the 4-body decay channel. 
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Figure 1: The ti becomes the NLSP in the whole marked region in the ~ tan/3 plane in the 
mSUGRA model for mo = 200(140), mi /2 = 145(180), // > 0 and mt =172.7 in the left (right) 
panel. The bounds on sparticle masses and other SUSY constraints, see the text, are also imposed. 
The CP-even Higgs boson mass contour for 111.4 is shown and the region above the contour is 
allowed (see the text for a brief review of the theoretical uncertainties in the calculated m/jo). The 
whole marked region is ruled out by the experimental bound m/jO ~ 114.4. Several contours of total 
4-body decay BR lighter top-squark are shown by different point styles, see the legends. The ‘-I-’ 
shows the regions where BR{ti bxi^i'e) — 20% for i=e and ji. The regions marked with heavy 
dotted square box are excluded by the CDF limit ( ~ 102 CeV for m~o ~ 51 CeV) obtained 

by assuming 100% BR of the loop decay. 
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Figure 2: The ti becomes the NLSP in the whole marked region in the Aq — tan/3 plane in the 
mSUGRA model for mo = 200(140), mi /2 = 145(180), // > 0 and mt =178.0 in the left (right) 
panel. The bounds on sparticle masses and other SUSY constraints, see the text, are also imposed. 
The CP-even Higgs boson mass contours for 111.4 (111.4 and 114.4) are plotted in the left (right) 
panel and the regions above the contours are allowed from the respective bounds. The whole marked 
region in the left (right) are ruled out from ru/jO ~ 114.4 (117.4) (see the text). The total 4-body 
decay BR’s contours of lighter top-squark are shown by different point styles, see the legends. The 
‘-I-’ shows the regions where BR{ti bxi^t^e) — 20% for 7=e and fi. The regions marked with 
heavy dotted square box can be excluded by the CDF limit (see Fig.l caption) . 


18 







&(GeV) 



A(GeV) A(GeV) 


Figure 3: The (mho) a function of for for mo = 140, mi /2 = 180, fi > 0 and tan /3=5,10,15 
in the left (right) panel in mSUGRA for mt=178.0. The ti-NLSP criterion is relaxed in this 
particular figure. 
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Figure 4: The ti becomes the NLSP in the whole marked region in the mo — mi /2 plane in the 
mSUGRA model for Aq = —630, ^ > 0, mt=178.0 and tan/3 =5(9) in the left (right) panel. 
The bounds on sparticle masses and other SUSY constraints, see the text, are also imposed. The 
CP-even Higgs boson mass contour for 111.4 (114.4) is plotted in the left (right) panel. The m/^o 
happens to be ~ 111.4 in the whole marked region in the right panel. The regions above the contour 
are allowed from the respective bound on rrij^o . The whole marked regions in the left (right) panel 
are ruled out from m/jO ~ 114.4 (117.4). The contours of total 4-body decay BR of the lighter 
top-squark are shown by different point styles, see the legends. The regions marked with heavy 
dotted square box can be excluded by the CDF limit (see Fig.l caption). 
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Figure 5: The mj^{mf^o) as a function of in the MSSM in the left (right) panel for mt=172.7. 
The MSSM model parameters are: M2=300, ^=300, Ma= 500, ^^=200, ylT-=300, m~ = 400 and 
'iTij= 250. The choices of tan/3 are shown in the legends. The ti-NLSP criterion is relaxed in this 
particular figure. 
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Figure 6: The ti is the NLSP in the whole marked region consistent with other sparticle mass 
bounds in the Xr — tan/3 plane in the MSSM model for M2=300, //=300, Ma= 500, Ah=200, 
y4T-=300, m~ = 400, mj= 250 and rrit =172.7(178.0) in the left (right) panel. The total 4-body 
decay BRs of the lighter top-squark are shown by the contours with different point styles. The ‘+’ 
shows the regions where BR{ti ~ 20% for i = e and /r. The Higgs boson mass contours 

for 114.4, 116.0, 118.0 and 119.0 (114.4, 118.0, 120.0 and 122.0 ) are also plotted in the left (right) 
panel. 
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Figure 7; The ti is the NLSP in the whole marked region in the — Ma plane in the MSSM 
model for M 2 =300, 300, m^400, m^250, Af,=200, ^,-=300, rrit =172.7 and tan/3=5(9) in 

the left (right) panel. The total 4-body decay BR of lighter top-squark are shown by different point 
styles. The ‘-I-’ shows the regions where BR{ti bxi^i'e) — 20% for i=e and /U. The Higgs boson 
mass contours for 111.4, 114.4, 115.0 and 115.7 (116.0, 117.0, 117.4, 117.8, 118.4 and 119.0 ) are 
also plotted in the left (right) panel. 
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